A thermodynamic force generated by chemical gradient and adsorption reaction 
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Biological units such as macromolecnles, organelles, and cells are directed to a proper location 
under gradients of relevant chemicals. By considering a macroscopic element that has binding sites 
for a chemical adsorption reaction to occur on its surface, we show the existence of a thermodynamic 
force that is generated by the gradient and exerted on the element. By assuming local equilibrium 
and adopting the grand potential from thermodynamics, we derive a formula for such a thermo- 
dynamic force, which depends on the chemical potential gradient and Langmuir isotherm. The 
conditions under which the formula can be applied are demonstrated to hold in intracellular reac- 
tions. The role of the force in the partitioning of bacterial chromosome/plasmid during cell division 
is discussed. 



Introduction Biological entities such as cells, or- 
ganelles, and macromolecules often move to a proper 
location under an external gradient of chemicals. In 
the chemotaxis of a cell, the process by which an or- 
ganism senses the presence of an external chemical and 
responds to it has been elucidated in depth. However, 
such coordinated motion is not limited to the organism 
level. In recent studies, it has been revealed that such di- 
rectional motion under chemical gradients/localizations 
plays an important role in organization at an intra- 
cellular level e^., microtubule guidance under the 
RanGTP gradient or Stathmin gradient 0], actin nu- 
clcation under the IcsA gradient on the outer membrane 
of the pathogen Shigella flexneri (actin comet) bac- 
terial chromosome/plasmid DNA partitioning by a Par 
system H 0, 0, i, 11 [13 [IH, Ell- and several trans- 
port processes by scaffold proteins. However, the gen- 
eral mechanism by which organelles or macromolecules 
demonstrate coordinated motion under an intracellular 
gradient/localization is yet to be determined. 

The pattern formation process of chemical concentra- 
tions under a non-equilibrium condition, first reported 
for a macroscopic morphogenesis jl4j . has recently been 
observed at an intracellular level as well, and its rel- 
evance to intracellular organization processes has been 
demonstrated in bacterial plasmid DNA partitioning Par 
system0[l3|, Min system for determining the plane of 
cell division in bacteria [l^, etc. They discussed how 
positional information given by the chemical concentra- 
tion gradient/localization is obtained and stored with 
the mechanism. However, there are few studies on how 
the coordinated motion or transport of organelles or 
macromolecules caused by chemical gradient /localization 
is sustained under non-equilibrium conditions. 

In the present paper, we present a physical mechanism 
of such directed motion by showing that under a chem- 
ical gradient corresponding to a reaction, a macroscopic 
element consisting of a number of reaction sites is gen- 
erally subject to a force. By representing this element 
as a scaffold for the adsorption of a chemical, we derive 
a general formula for a force generated by a chemical 
potential gradient. In the derivation, we introduce the 



grand potential in thermodynamics for an open system 
jl7j and apply the second law of thermodynamics. Irre- 
spective of the nature of the molecules, the force acts to 
decrease the grand potential when a chemical potential 
gradient exists. Then, the element moves in a direction 
such that the chemical potential increases. A formula is 
obtained by assuming that the reaction process reaches 
equilibrium faster than the motion of the element and 
by extending the minimization of free energy to include 
the contact with a particle bath with a given chemical 
gradient. We propose that this force generates a gen- 
eral mechano-chemical coupling driven by the chemical 
potential gradient through an adsorption reaction. The 
proposed formula is valid since the work extracted by 
the gradient force dominates the thermal energy fluctua- 
tion. By examining whether this condition is satisfied in 
the intracellular reaction process, we discuss the possible 
role of this force in the partitioning dynamics of bacterial 
chromosome/plasmid for cell division. 

Chemical gradient force Now, consider organelles or 
macromolecules that have a number of binding sites for 
reactions to occur on their surface; for example, nucleo- 
protcin complexes (NCs) have several promoter sites to 
bind transcription factors, etc. Let us model these bi- 
ological elements simply as "beads" with several reac- 
tion sites to which the corresponding molecules attach, 
as shown in FiglTJ The bead is placed at r = ^ and 
moves in a d-dimensional space r G (d=l,2,3). We 
consider an isothermal process that is homogeneous over 
space at a given temperature T , whereas we consider a 
chemical bath with a spatially dependent concentration 
x{r) of a chemical X, or equivalently, the corresponding 
chemical potential fi{r). This gradient is assumed to be 
sustained externally. This chemical is attached to the 
binding sites B on the bead and thus forms a complex 
Y (Figd]), as given by the reaction X(^) B ^ Y. The 
molecular number of the complexes on the bead is de- 
noted by Ny. Note that the "bead" here is defined as a 
macroscopic entity compared to the molecule X. To de- 
velop local-equilibrium conditions, we make the following 
assumptions. The "adsorption" reaction on the bead is 
considered to be a macroscopic event and Ny is a molec- 
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ular number averaged over a much longer timescale than 
the microscopic timescale of the reaction. Further, the 
bead is assumed to move sufficiently slowly so that lo- 
cal chemical equilibrium of the above reaction can be as- 
sumed to exist at each position r ~ ^. In other words, the 
timescales of diffusion of X molecules Tdi / / and adsorp- 
tion reaction Tadsorb are much smaller than the timescale 
of the motion of the bead n^ad- Tdif f ,Tadsorb < Uead- 
As long as the bead is in motion, it is maintained in 
equilibrium with the reservoir at r = ^. 

With the assumption of the existence of local equilib- 
rium, we can apply thermodynamics with spatially de- 
pendent thermodynamic variables. Indeed, at each po- 
sition ^, the reaction process is described by the famil- 
iar classical Langmuir adsorption theory Note 
that the theory has been successfuly applied to DNA- 
protein binding equilibrium [l5| . The grand potential in 
thermodynamics at each position is given by ^{$,) = 
F{i) ~ vKil d^) = dF{i) ~ d{y,,m = -ydfiiO, 
where y = ^ {V : volume of the bead), F{$) is the 
Hermholtz free energy, and dF{$,) = f-i{S,)dy |17| [22|. 
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Figure 1: Schematic representation of our system. A "bead" 
is placed at r = ^ and moves in a d-dimensional space 
(d=l,2,3). The adsorption reaction X(^) + B Y occurs 
on the surface of the bead. 

Now, consider a virtual displacement for the position 
of the bead r = ^. Under an infinitesimal displacement 
^ ^ + d^, the change in the grand potential is = 
—ydfi^^) = — yV/i(^) -d^. In other words, the position of 
the bead ^ is adopted as an effective independent variable 
instead of the chemical potential ^(^). Then, ^ is the 
work coordinate, while — ?/Vju(^) is the force exerted on 
the bead by the external world balanced by the force 
exerted on the bead due to reservoir chemical potential 
distribution /i(r). The chemical gradient force generated 
by the reservoir acts on the bead per unit volume of the 
bead: 
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(1) 



This expression is interpreted as follows: consider a 
quasi-static infinitesimal displacement d^. Then, from 
the change in the grand potential, the maximum work 
that the system docs on the external world through the 



reservoir is d'W = yV/i(^) • d^. By considering that this 
work is due to the force that the reservoir exerts on the 
bead, as d'W = fchem ' the force formula cq.(l) is 
obtained. In other words, without an externally applied 
force, ^ evolves spontaneously so that monotoni- 
cally decreases, that is, dJl(^) < [2^. When we con- 
sider the overdamped system, where the kinetic energy 
of the bead is negligible, the equation of motion is given 
by = Tcr{^) > 0, where cr(|) := ^(£) is the en- 

tropy production of the system. The above equation of 
motion is rewritten as 7 ^ = — Vfi(^) = when 
there is dissipation due to friction, assumed to be propor- 
tional to the velocity with the proportionality constant 
7 — the friction coefficient per unit volume of the bead. 

Now, consider the condition for local chemical equilib- 
rium. Because c ~ y + b = const., the dissociative con- 
stant is defined as A' := = ElUt^ where y is given by 



the Langmuir isotherm y = yix{^)) — c-^j^^^ [2 
we consider cooperative adsorption of n chemicals on the 
bead: nX(^) -I- B ^ Y, the isotherm for the adsorption is 
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where n determines a Hill 



given as y{x{^)) = c-^„+^(^)„ ; 
coefficient. Therefore, the equation of motion is written 
as 
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(2) 



Equation (2) is a general formula for the motion of an 
element that has a number of binding sites for chemical 
adsorption under the gradient of corresponding chemical 
potential. The theoretical description of the motion of 
the bead is somewhat an extended form of Langmuir ad- 
sorption theory. Because of this force, the bead moves 
in a direction such that the grand potential f2, the ther- 
modynamic potential for open systems, is decreased; this 
decrease follows from the second law of thermodynamics. 
Accordingly, the bead moves in the direction of increas- 
ing chemical potential. 

In our formulation, the chemical potential fJ-{r) is re- 
garded to be equivalent to the external potential field 
for bead motion, similar to the magnetic field in which 
a magnetic dipolc is placed. Indeed, the magnetic force 
that acts on the dipolc placed in the magnetic field gra- 
dient acts in the direction of increasing magnetic field 
strength. 

When an additional external potential field E{r) such 
as elastic energy is applied to the bead, the relation 
dF{^) = dE{^) + n{k)dy{i) holds from the first law 
of thermodynamics. Accordingly, the grand potential 
of the system given by il(^) = F{^) — uiOf^i^) satis- 
fies dn{^) ^ dFi^) - d {Ki)ym = dE{^) - yiOd^- 
Further, by taking into account thermal fiuctuation, the 
equation of motion of the system is given by 



7^ = -Vi?(0+2/(OVM(|)+r7(t) 



(3) 



with (T7(t)) = 0, {r]{t) ■ r]{t')) = 2d'ykBTS{t - t'). Note 
that the equation is obtained as an extended form of 
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a conventional one with the Hermholtz free energy as 
the thermodynamic potential for a closed system. It is 
straightforward to further extend the present formula to 
include a case involving multiple beads and multiple com- 
ponents with interactions among the beads. 

An example We consider a simple toy model as 
an example of eq.(3). We assume that a bead placed 
at r = ^ is balanced at r = in a 1-dimensional 
space r by a restoration force, which is produced by a 
linear spring and represented by a harmonic potential 
-E(^) = f^^. The concentration of the chemical that 
reacts with the bead is assumed to be sustained and is 
given by x{r) = Xg exp(A(r — r^)), where Xs =x{rs) is the 
concentration of X at r = rg, while the corresponding 
chemical potential is given by /i(r) ~ Jl + kBThix{r). 
Here, p, is the standard chemical potential. Then, by us- 
ing eq.(3) without considering thermal fluctuation, the 
equation of motion of the bead is obtained as 7 C = 
^(0^(0-^(0 = AfcBrc ^„f4„ -ag. Here,7,c,iC,a 
are the frictional coefficient, maximum adsorption con- 
centration, dissociative constant, and spring constant, 
respectively. This equation can be rewritten as 7 ^ = 

-ttr T^-rn ^ - < with A = XkgTc^B = ( —Y ■ 

The solution of the equation shows bistability as it in- 
cludes a sigmoid function. Thus, under appropriate pa- 
rameter regions, as Xs, which is taken as a control param- 
eter, is increased, the model thus undergoes two-fold bi- 
furcations from monostable fixed points to bistable ones 
and then back to monostable states. This implies that 
the gradient force determines the stable fixed position of 
the bead. This toy model can be used to explain how the 
stable position of an organella undergoes bifurcation to a 
new location as the chemical concentration increases. An 
example may be seen in bacterial chromosome/plasmid 
partitioning, as will be discussed later. 

Conditions for chemical gradient force (i) : relation- 
ships among timescales The grand canonical descrip- 
tion with local chemical equilibrium can be adopted only 
in certain conditions. The following relationships among 
timescales must exist: (1) the timescale of adsorption 
Tadsorb must bc shortcr than that of bead motion Thead 
in order for the adsorption reaction to reach equilib- 
rium (approximately) before the bead moves, and (2) the 
timescale of bead motion must be much larger than that 
of the diffusion of the chemical adsorbed on the bead, 
i.e., Tbead 3> Tdiff- This condition is necessary to assume 
macroscopic motion of the bead under the chemical gra- 
dient, which itself diffuses in space. In the living cells of 
particular interest, the molecular number is not so large 
that the instantaneous chemical concentration fluctuates 
greatly. Hence, in order to determine the macroscopic 
concentration, the temporal average of the concentration 
of molecules over the diffusion timescale must be calcu- 
lated p^ . and hence, the above condition is important. 

Condition (ii): dominance of chemical gradient force 
over thermal fluctuation As long as the above con- 
ditions are satisfled, on an average, the gradient force 



acts in the direction in which the chemical potential 
increases, regardless of thermal fluctuation. However, 
the gradient force must be larger than thermal noise in 
order for the force to act effectively without the need 
for long-time averaging to remove fluctuation. This 
establishes the following condition: the work done by 
the force must be larger than the thermal energy, i.e., 
—V c?C'2/(€')^m(€') ^ kgT for directional motion 
from £q to ^ [2^ . In the case of one-dimensional motion 
from 7' = to L, the above condition for the work is 
rewritten as follows: the concentration and chemical po- 
tential of chemical X arc obtained as x{r) = xq exp(— ^) 
and fi{r) ~ p, + kBTliix{r), respectively. Here, xo 
gives the highest concentration at r — 0. Then, the 

above inequality is written as /q 0?^ — '—^ > 1 

with Nc = Vc. A straightforward calculation leads 



to the inequality — In 



1+ 



(^f)"exp(-ni) 



> 1. 



can rewrite the above inequality in the form — 
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with the additional condition 



Note that the the lower bound on ^ 
and the formula is valid over 



1— oxp( — n( 

decreases with — and N, 
a wide range of ^ . 

An application Now, we consider applications of the 
chemical gradient force within a cell. As an example, 
consider the partitioning of chromoso me /p lasmid in bac- 
teria during cell division [1, % 0, i, i, [13 [O, [11 • Here, 
the bead corresponds to the chromosome/plasmid on 
which the relevant protein (X in the model) binds to 
form a nucleoprotein complex (NC) (Y in the model) 
that is important for partitioning. For the applica- 
tion, the above-mentioned condition on the timescales, 
Tdiff, Tadsorb < Uead has to bc satisflcd. In general, 
the timescale of protein binding equilibrium on bacte- 
rial DNA Tadsorb IS at most Tadsorb ^ l(s)- In , it 

has been suggested that the segregation of chromosomes 
is spatially restricted and the diffusion coefficient D is 
estimated as O {iQ"^ - 10"^). The size of NC, a, is as- 
sumed to be on the order of a 50(nm), so that Th^ad 
is estimated as Uead ~ 25 (s). Similarly, 

Tdiff is estimated as Tdiff ~ 10 '^(■s), where is 

the diffusion coefficient of proteins within cytoplasm and 
is roughly estimated as ~ 3.0(^m^/s) (according to 
[131). Therefore, Tdiff < Tadsorb < Uead is satisfied for 
our formula. 

To be specific, the replicated bacterial chromo- 
somes/plasmids are activelypartitioned to daughter cells 
by a Par system 0, i, 0,lS B [13, El, [ll- The most 
ubiquitous partitioning system is the parABS system. 
Before cell division, daughter chromosomes/plasmids are 
precisely segregated by ParA dynamic localization along 
with the associated ParA-NC interaction. It has been 
suggested that ParA localization pulls chromosomes or 
at least guides plasmid [l^l, but how ParA generates the 
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driving force for chromosomc/plasmid segregation has 
not been elucidated yet. Now, ParA is regarded as X, NC 
(to be precise, ParB-par^* NC mediated by ParA ) as Y 
in our model [23] ; this driving force is naturally explained 
by the chemical gradient force. To confirm the validity 
of our formula, we simply consider the steady concen- 
tration distribution of ParA to be x{r) = a;oexp(— ^) 
within [0,L] 0,0]. We then take n = 2, because ParA 
cooperatively binds DNA [2l| and choose K = 0.3(^M) 
on the basis of recent data obtained by [2l| . Although the 
precise number of ParA that binds DNA around NC is 
unknown, Nc 10 is a natural choice since many ParA 
are suggested to bind DNA around NC. The localiza- 
tion scale Tc of ParA localization is also not precisely 
known and could vary across bacterial species. By tak- 
ing ~ 5 as an order estimate from [g, Q, we obtain 
^ > 0.5. Active protein binding generally occurs when a 
reasonable concentration xq K is reached, where ^ is 
sufficiently large to satisfy the above inequality. There- 
fore, the chemical gradient force dominates the thermal 
fluctuation if xq > K. Although the mechanism of ParA 
localization could depend on the species, the estimate 
suggests that the force plays a significant role in the 
partitioning of plasmid/chromosome, independent of the 
specific molecular mechanisms. Further, wc note that 
the above-mentioned "spring toy model" under a gradi- 
ent is applicable to the partitioning problem and may 
help explain the bifurcation of plasmid/chromosome in 
cell division with the change in ParA localization. 

Of course, the biological ParA distribution is more 
complicated. Although the application mentioned here 
assumes ParA steady distribution, the distribution some- 
times shows pole-to-pole oscillation [l2| and regular pat- 



terns to form regular spacing of a few plasmid foci 
[13, Eli- Still, our theory is applicable to the time- varying 
distribution, because the force derived here leads to di- 
rectional motion toward the ParA foci. Indeed, the di- 
rectional motion toward ParA focus has been suggested 
in [l^ . Furthermore, the present gradient force is appli- 
cable to reaction-diffusion system so that the plasmid is 
located according to the chemical concentration pattern. 

Discussion In this study, we have demonstrated that 
when a chemical concentration gradient exists, a macro- 
scopic element that acts as a scaffold for the absorption of 
some chemical is generally subject to a thermodynamic 
force, and it moves in the direction of increasing chemi- 
cal potential. The force is generated in accordance with 
the second law of thermodynamics, as represented by the 
decrease in the grand potential in thermodynamics un- 
der a given gradient; thus, this force acts independent 
of the specific molecular mechanism. The gradient force 
is expected to provide a novel perspective on mechano- 
chcmical coupling. We have derived a general formula 
for the magnitude of the force and presented a few con- 
ditions necessary for this formula to be valid. The condi- 
tions are shown to be satisfied for the intracellular motion 
of macromolecules or organelles under a suitable chemi- 
cal gradient, and the force is clearly shown to be larger 
than the thermal fluctuation force. To be specific, we 
have discussed the application of this formula to bacte- 
rial partitioning systems in cell division. Considering the 
generality of our formulation, the chemical gradient force 
is expected to be a universal guiding principle in intra- 
cellular organization. 

The authors would like to thank H.Takagi, A.Awazu, 
S.Ishihara, and T.Yomo for discussions. 
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proved by the principle of maximum work in an isother- 
mal process. See Supplementary Information for a de- 
tailed derivation. 

[24] In terms of statistical mechanics, the Langmuir isotherm 
should be obtained as a function of the chemical potential 
l-L{r) at r = ^. We can derive the representation as a 
function of ^i(r) when ^t(r) is that for a dilute solution, 
i.e., /i(r) = fi + kBTlnx{r). Here, fi is the standard 
chemical potential. 

[25] This inequality is also derived from the steady solu- 
tion of the Smoluchowski equation. For one-dimensional 
motion, the steady-state distribution of ^ is given by 

P{£,) = Aexp ^— ^ ■ Here, A is a normalized fac- 
tor and W{0 = -V"/(f dC'y(5')VAi(?')- ^ the position 
that satisfies W{(,c) ~ ksT is smaller than the sys- 
tem size L, the distribution is concentrated toward higher 
chemical potentials and the gradient force dominates the 
thermal fluctuation. Then, the inequality in the text is 
obtained. 

[26] We can also represent the condition in terms of the av- 
erage concentration Xav for [0,L]. 

[27] The Par system consists of three components: DNA 
binding protein ParB, ATPase Par A, and centromere- 
like site parS. ParB binds parS, spreads along DNA, 
and forms a large NC around parS. ParA can nonspecif- 
ically bind to DNA and interact with the ParB-parS NC 
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Supplementary information 

Proof of the variational principle for grand potential 
with the principle of maximum work — It is well known 
that the grand potential is a thermodynamic potential in 
terms of natural variables T, and surface area; exten- 
sive studies such as those on adsorption and capillary 
phenomenon involve this potential. However, to the best 
of our knowledge, there is no proof of the variational prin- 
ciple for the grand potential, particularly when fi varies; 
nevertheless, this potential is closely connected by the 
Legendre transformation to the other thermodynamic po- 
tentials for which the existence of variational principles 
has been proven. It seems to be because 

• standard thermodynamic operation involves exten- 
sive variables but not intensive ones in order to ob- 
tain work from a thermodynamic system. 

This is not the case in this study because fi can be con- 
trolled by changing the d-dimcnsional position of the 
thermodynamic system {d— 1,2,3). Therefore, we will de- 
rive the variational principle for the grand potential on 
the basis of the principle of maximum work. 

Consider a thermodynamic system composed of a bead 
with a surface as shown in (Fig.l). The bead is placed 
at r = ^ and moves in a d-dimcnsional space r G R 
{d— 1,2,3). It makes contact with a heat-particle bath 
(reservoir) with homogeneous temperature T{r) = T and 
chemical potential /z(r) (concentration x{r)) at r ^ 
for a chemical X. The chemical is adsorbed on the sur- 
face of the bead B and forms a complex Y with it 
(Fig.l): X(^) -I- B ^ Y. The concentration on the bead 
is denoted by y. 

Assume that the bead is in equilibrium with the reser- 
voir at r = £ and it is maintained in equilibrium during 
its motion. Therefore, the temperature and chemical po- 
tential of the system are always balanced by those of the 
reservoir, that is, Tgystem = T,iisystem = Ki)- Note 
that we assume the existence of local equilibrium, where 
several thermodynamic variables can be defined in the 
d-dimensional space. 

The precise definition of the grand potential is as fol- 
lows: 

^{T, tj) ■■= min [F{T, y) - yfi] 
y 

This is the precise representation for the Legendre trans- 
formation from F{T,y), the Helmholtz free energy in 
terms of natural variables (T,y), to Q{T,fi) in terms 
of natural variables {T,ii); mathematically, this implies 
that y must be determined as a function of T, fi, so that 
F(T, y) — yfi is minimized when y is varied with fixed 
T, /i. The determined y* and F{T, y*) — y*^ are iden- 
tical to the Langmuir isotherm and the grand potential, 
respectively. 

The principle of maximum work represents the second 
law of thermodynamics in terms of the work obtained by 



an isothermal process. First, we assume that the princi- 
ple of maximum work is applicable, and we vary y while 
keeping T constant: 

WmaAT, yi ^ y2) = F{T, y,) - F{T, y^) 

When y is quasistatically changed from j/i as a func- 
tion of /ii to any y2 by moving the position of the bead 
from to ^2 via the force / exerted by the external 
world (This force was balanced by the chemical gradient 
force or forces from the other potential fields, such as me- 
chanical or magnetic forces.), the work W{T,^i ^2) 
in terms of natural variables T, fi obtained by changing 
fi from fii to fi2 and keeping T fixed, which is the work 
performed by the system on the external world and is 
given as follows: 

hi 

= Wynax{T, 2/1 2/2) - f - / d (yn) j 
V Jyil-i-i ' 

The first term represents the total work that the sys- 
tem can perform, and the second term represents the 
work that it must perform on the reservoir for the change 
2/1 ~* 2/2, i-c, the work that it cannot perform on the ex- 
ternal world. We can interpret the latter work as follows: 
the potential energy of the system is — 2/(^)/^(^) at r = ^ 
when it is placed in a potential field /i(T"), similar to a 
situation in which a dipole is placed in the magnetic field 
generated by a magnet, as discussed below. Furthermore, 
the former (latter) work is related to the interactions that 
result in heat (chemical X) exchange with the reservoir, 
respectively The interactions are represented by chang- 
ing y through ^ while T remains fixed. The essential 
point is that only the latter work yields information on 
the system as an open system. 

We consider maximizing W{T,yL\ fi-i) by changing 
2/2 determined arbitrarily so far and then determining it 
as a function of T, [12. The maximum work is expressed in 
terms of natural variables T, fj, as Wmax{T, fii -^^12). We 
obtain the maximum work in terms of natural variables 
(T, (i) as follows: 

Wmax{T,ijLi -> ^2) = maxH/(T,yLti ^ 112) 
y2 

= max [WmaxiT, 2/1 2/2) + 2/2/^2 - 2/iMi] 
y2 

= F{T, /ii) - 2/iA*i - min [F{T, ^2) - 2/2^*2] 

= ^2(r,Mi)-o(r,/i2) 

That is, the maximum work that the system assigned 
(T, /i(^)) by the reservoir {T,fj,{r)) at r = ^ can per- 
form on the external world is obtained by calculating the 
decrease in the grand potential fl{T, /i($)). 
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Next, we consider the situation in which the force / 
exerted by the external world vanishes (/ = 0) and 
the change /i — > occurs spontaneously. In this case, 
W{T, ^')^0 because / = 0. 

.-. w^ma.(r,// ^ fi') > M^(r,/x -> = o 
n{T,fi') - n{T,fi) < 

Because the start state (T, fi) can be chosen arbitrarily, 
the grand potential monotonically decreases till the end 
state (T, /x') under local equilibrium is reached. If (T, fi') 
is in equilibrium, the grand potential reaches a minimum 
at this state. 

For an infinitesimal displacement ^ ^ ^ + d^, the fol- 
lowing inequality is satisfied: 

dn{T,ti) = vo(r,o -rfC < 

Therefore, for a displacement per unit time ^, 

= vo(T,o-^<o 

Because the change in the grand potential is completely 
dissipated due to W{T, fi ^ fj, + dfi) = if the kinetic 
energy of the bead is negligible, the following inequality 
is satisfied: 

-f2(r,ei) = ra(i) >o 

ct(^) is the entropy production of the system. Thus, we 
can show that the variational principle applies to the 
grand potential of the system. 

If the dissipation in the system is only caused by fric- 
tion due to the motion of the bead, the following relation 
is satisfied for the power that the system does for friction 

-n{T,t^)^Ta{^)^Wfri^)>0 



Comparison with Gibbs free energy of magnetic dipole 
moment — For understanding the work that the sys- 
tem must perform on the reservoir as discussed above, it 
may be convenient to compare the problem with that of 
a magnetic dipole moment placed in a magnetic field. 

It is well known that a force exerted on the dipole 
placed in a magnetic field acts in the direction of the 
larger one when the strength of the magnetic field is in- 
homogeneous. The Gibbs free energy of the magnetic 
dipole moment in terms of natural variables T, H, for 
which its demagnetizing field is negligible, is defined as 
follows: 

G(T, H) := mm[i^(r, M) - MH] 

F{T,M),H,M are the Hclmholtz free energy of the 
dipole with natural variables T, M , magnetic field, and 
magnetization per unit volume, respectively. With T 
fixed, 

dG^dF - d{MH) = -MdH = -MVH ■ d| 
Now we will show that dG{T, H) < 0, similar to the 
above discussion on that of the grand potential. When 
the dipole at r = ^ is placed in the magnetic field H{r), 
it has potential energy —M{^)H{^). This means that 
under an infinitesimal displacement ^ ^ + d^, the sys- 
tem must perform work —d {M{^)H{^)) on the magnetic 
field in addition to the work performed on the external 
world. It is clear that this is equivalent to —d{y{^)fj,{^)), 
when y, fi are replaced with M, H in the above discussion. 
Therefore, we see that when T is fixed, dG{T, H) < 0. 
When there is dissipation due to friction, assumed to be 
proportional to the velocity with the proportionality con- 
stant 7, the equation of motion for the dipole is given as 
follows: 

7^ = MVH 



